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1. Motivation 

Twisted mass lattice QCD (tmLQCD) offers an efficient mechanism for eliminating unphysical 
zero modes[[T]] and removing 0{a) en-ors[2] from Wilson simulations. Given the successes of 
tmLQCD in dealing with the {u,d) quark doublet, it now remains to extend the theory to include 
other quark flavours, noting that the next lightest pair, {c,s), is far from being degenerate. The 
present work makes a step in this direction by investigating one option. 

Some unphysical flavour symmetry violation occurs at non-zero lattice spacings in tmLQCD. 
For example, in the two flavour theory with a degenerate (i<, J) doublet, [||, ^] 

(niA++ = mA- ) / {m^o = ) . (1.1) 

When tmLQCD is extended to include the strange quark, many splittings arise among strange 
hadrons as well. Since all unphysical flavour splittings must vanish in the continuum limit, they 
provide a measure of discretization enors in tmLQCD. 



2. Lattice action 

There is more than one way to incorporate the strange quark into tmLQCD. In the present 
work, we use a tmLQCD fermion action similar- to that of Ref. [§] which is a four-quark action, 
though the fourth (charm) quark will not affect any of our quenched three-flavour spectrum com- 
putations. The action is 

Sf[Y,V,U]=a^Lvi^) (^ + 'At75 + 7-V^-^£v;Vv) vi^), (2.1) 

where 

= {u,d,c,s), (2.2) 

M = diag{Mi,Mi,M„M,), (2.3) 

H = diag{ni,-Hi,Hc,-l^s)- (2.4) 

The parameters M„ and n„ denote the standard quark masses and twisted quark masses respectively. 
Forward (V), backward (V*) and symmetric (V*) derivatives are standard. To obtain automatic 
0{a) improvement we choose maximal twist. The classical definition would be 

(Oi = arctan^^^ = |, (2.5) 

«2^arctan(^)=arctan(^)=|, (2.6) 

but quantum corrections lead to renormalization of both M and Maximal twist is therefore 
defined nonperturbatively by||3|, |^, ^ 

^ ^X.v(Vv-(x,0P+(0)) ^ 

tana)„ = jhz — - — ; — (2.7) 

where a tilde denotes currents constructed from fields in the twisted basis. The correlators on the 



right hand side of Eq. (2.7) require only a single fermion propagator, hence the single subscript n. 



In the present work, this fermion action will be combined with the standard Wilson gauge action. 
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Table 1: Parameter choices for the simulations. 

3. Simulation details 

Table [T] shows the parameter values chosen for our action and lattices. The authors of Ref. [^] 
used ro to determine the lattice spacings: a = 0. 123 fm for j8 = 5.85 and a = 0.093 fm for j8 = 6.0. 
The standard mass parameter values, aMq, were determined in Ref. [||] to correspond to maxi- 



mal twist as defined by Eq. (2.7). The twisted mass parameters range from the approximate physical 
strange quark mass, nis, to 1/6 of m^. Notice that the standard mass parameter is tuned for each 
individual twisted mass parameter so, for example, the quark and anti-quark in a computation of 
the kaon mass will have separate aMq values. 

Local interpolating operators ai^e used throughout this work, and three state fits to the resulting 
correlators are performed. 



4. Results for meson masses 



Numerical results for ground state pseudoscalar meson and vector meson masses are displayed 
for j3 = 6.0 in Fig. |l]. (Comparable results were obtained for j8 = 5.85.) All quark/anti-quark mass 
combinations from Table |I] are shown in this figure, with an^ows indicating the values of twisted 
quai^k masses (i.e. the horizontal axis of each plot) corresponding to the physical hadrons. For 
pseudoscalar mesons, significant mass splittings are evident between charged and neutral mesons, 
as read from the vertical axis of each plot. 

The observed mass splitting between TT^ and 7r° is affected both by quenching and by our 
omission of "disconnected" contributions. (See Ref. |^] for a recent tmLQCD study of this split- 
ting.) However, the to K^^ splitting requires no disconnected contributions and thus our 
data provide the complete quenched result in that case. Given the automatic 0{a) improvement 
of masses in tmLQCD ||^, this unphysical kaon splitting should vanish linearly with a^. Figure ^ 
shows that our two j8 values are consistent with this expectation. 

In contrast to the pseudoscalar mesons. Fig. |l| shows no significant splitting between charged 
and neutral K* mesons within the statistical uncertainties. 



5. Results for baryon masses 

The lack of flavour symmetry in tmLQCD (with a^O) means the physical octet of spin- 1/2 
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Figure 1: Pseudoscalar (left panel) and vector (right panel) meson masses as functions of the sum of quark 
and anti-quark twisted masses at j3 = 6.0. 
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Figure 2: Mass splitting between neutral and charged kaons as a function of squared lattice spacing. 



baryons has no exact mass degeneracies and the decuplet of spin-3/2 baryons has just two:[§, 

mA++ = m^-, (5.1) 
ni^o = (5.2) 

For our numerical study of baryon masses, the local source and sink operators are given in the 
legends of each panel in Fig. ^. All combinations of quark masses from Table [I] were used at each 
j3 value for both spin- 1/2 and spin-3/2 operators, though Fig. ^ only displays data from one j3 as 
an example for each spin. 

Any flavour splittings among these numerical results for baryons are at the edge of statistical 
significance. Clarification of these splittings will require better operators and increased statistics. 

Figure]^ compares the baryon masses in physical units to the spectrum observed in nature. The 
tmLQCD results lie systematically above experiment, which could be a quenching artifact and/or 
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Figure 3: Panels (a), (b) and (c): (1/2)+ octet baryon masses, plus one (1/2)^ channel to discuss the 
A(1405), as functions of the sum of the twisted quark masses at j8 = 6.0. Panels (d), (e) and (f): (3/2)+ 
decuplet baryon masses as functions of the sum of the twisted quark masses at j8 =5.85. 



a consequence of extracting the lattice spacing from tq. The two jS values produce essentially 
equivalent results, suggesting that we are already in the scaling region. 



6. Conclusions 

One particular method for incorporating the strange quark into tmLQCD has been explored 
by directly computing the hadron mass spectrum. This approach to 3-flavour tmLQCD produces 
a large splitting among kaon masses at a > 0.1 fm. The splitting vanishes linearly with as 
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Figure 4: The quenched baryon spectrum compared to experiment. 

c? ^ as expected. Some of the anticipated baryon splittings are potentially discernable, though 
small. Because the fermion determinant is complex, this method is better suited for quenched 
simulations than for dynamical ones. Other options for tmLQCD with a strange quark should also 
be investigated and compai^ed. 
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